A maximally superintegrable system on an 
n-dimensional space of nonconstant curvature 



A. Ballesteros a * A. Enciso 6t F. J. Herranz * O. Ragnisco c§ 

a Depto. de Ffsica, Universidad de Burgos, 09001 Burgos, Spain 

b Depto. de Ffsica Teorica II, Universidad Complutense, 28040 Madrid, Spain 

Dip. di Fisica, Universita di Roma 3, and Instituto Nazionale di Fisica Nucleare 
Via Vasca Navale 84, 00146 Rome, Italy 



Abstract 

A novel Hamiltonian system in n dimensions which admits the max- 
imal number 2n — 1 of functionally independent, quadratic first inte- 
grals is presented. This system turns out to be the first example of a 
maximally superintegrable Hamiltonian on an n-dimensional Riemannian 
space of nonconstant curvature, and it can be interpreted as the intrinsic 
Smorodinsky-Winternitz system on such a space. Moreover, we provide 
three different complete sets of integrals in involution and solve the equa- 
tions of motion in closed form. 



1 Introduction 

A Hamiltonian system on a 2n-dimensional phase space is maximally super- 
integrable (MS) if it is Liouville integrable [2] and there exist further n — 1 
functionally independent (global) first integrals. Throughout this paper we 
shall restrict ourselves to the case of natural Hamiltonians, whose kinetic term 
defines a Riemannian metric on an n-dimensional configuration space. 

The intrinsic relevance of MS systems has resulted in a wealth of results con- 
cerning MS Hamiltonians [11,20,21,28,34]. MS models with integrals of motion 
quadratic in the momenta are of particular interest because of their connections 
with generalized symmetries [19,35], isochronic potentials [10,18], and separa- 
bility of the associated Hamilton-Jacobi (HJ) equation and Schrodinger equa- 
tions [8,22,27]. A systematic study of MS systems on certain low-dimensional 
spaces is currently being developed by Kalnins et al. (cf. [24] and references 
therein), and in fact in two dimensions they managed to obtain a (local) clas- 
sification of all MS Hamiltonians with integrals at most quadratic in the mo- 
menta [25,26]. This result, which hinges on a classical theorem by Koenigs [29], 
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yielded the first examples of MS systems not defined on constant curvature 
spaces. 

Unfortunately, most examples of MS Hamiltonians are restricted to low- 
dimensional spaces. To the best of our knowledge, the only known examples of n- 
dimensional MS systems with quadratic integrals are given by the (generalized) 
Kepler problem [14], the Smorodinsky-Winternitz (SW) system [15,17], their 
generalizations to the simply connected spaces of constant sectional curvature 
(cf. [4] and references therein), and the geodesic flow on these spaces. Note that 
the allowed potentials correspond essentially to those appearing in the extension 
to space forms of Bcrtrand's theorem [7,31]. If one drops the requirement 
that the integrals be quadratic in momenta, the list can be enlarged with the 
rational and hyperbolic Calogero-Sutherland-Moser models of type A [9,32,36, 
38], the nonisotropic oscillator with rational frequencies and the nonperiodic 
Toda lattice [1,37]. 

In this paper we show that the Hamiltonian 

H(p ' q) - ^ ' (1) 

where n > and uj 2 , b } > 0, is a MS system with quadratic first integrals. Here 
p, q G W 1 , and we shall always assume that R 2n is endowed with its standard 
symplectic form dq A dp. For any u £ M. n we use the notation 

u = u • u = Ui 

i 

and all the sum indices run from 1 to n unless otherwise stated. 
The Hamiltonian {!]) yields the motion under the potential 

V q = /V J ( 2 ) 

k -\- q 

on the conformally flat Riemannian manifold M. n := (R™,ds 2 ) whose metric is 
given by 

ds 2 = (K + q 2 )dq 2 . (3) 

The Hamiltonian |T]) provides the first example of a Hamiltonian system on a 
Riemannian space of nonconstant curvature which is MS in any dimension. The 
analysis of the geometry of M, n presented in Sec. Q] shows that this model can 
be in fact regarded as the intrinsic SW system on this space. When n = 2 this 
model is listed in Kalnins et al.'s classification of MS systems in the Darboux 
space of type III [26] . This model is singled out among the others in this list 
because it is the only one readily amenable to the Poisson coalgebra treatment 
developed in Sec. O 

The paper is organized as follows. In Sec.[Uwe show that the Hamiltonian ([T]) 
is indeed MS and briefly study the separability of the HJ equation. The proof 
is greatly simplified by the use of the fit(2) Poisson coalgebra symmetry [5,6] of 
the model. In Sec. [3] we study the equations of motion and find their general 
solution in closed form using the previously computed first integrals. Quite 
remarkably, the radial motion is essentially that of the Kepler problem, whereas 
the HJ equation essentially coincides with the one for the SW system. In Sec. |4] 
we analyze the geometry of the underlying Riemannian space and relate the 
Hamiltonian fTj to the intrinsic harmonic oscillator in A4 n . 
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2 Maximal superintegrability and separability 
of the HJ equation 

Instead of dealing with the model (p} directly, in the following two sections we 
prefer to analyze the Hamiltonian 



p 2 -c + Ejfe jgr 2 

2{k + q 2 



H(p,q) = - nf ,\C 3 , c:=klo\ (4) 



which is related to ([TJ) via Ti — 2H + lj 2 . Obviously the maximal superinte- 
grability of H is equivalent to that of H, whereas the trajectories of H can be 
mapped into those of H via the time dilation t \— ► t/2. In particular, H and TL 
possess the same orbits. 

A first observation is that there is an st(2) Poisson coalgebra symmetry of 
the Hamiltonian ((4| that immediately provides 2n — 2 quadratic integrals. This 
underlying coalgebra structure is shared by all the other known examples of 
quadratically MS models [3,4], as well as by uncountably many other integrable 
Hamiltonian systems which are not MS. Therefore, the crucial result in this 
section is the derivation of an additional first integral, functionally independent 
of the previous 2n — 2. 

The s[(2) Poisson coalgebra is defined by the basis { J £ : e = ±, 0} together 
with the following Lie-Poisson brackets, Casimir C and primitive coproduct A: 

{J ,J+} = 2J+, {J , J_} = -2J_, {J_,J+} = 4J , (5) 
C = J-J+- J 2 , A(J e ) = J e ®l + l®J £ . (6) 

A symplectic (one-dimensional) realization of the algebra ([5|) on (M 2 , d<?i A dpi) 
is given by 

j- = Qi, Jo = npi , j+ = p\ + Mr 2 . 

where b\ is a real parameter labeling the representation which corresponds to the 
value of the (one-dimensional) Casimir. Under this realization, the Lie-Poisson 
product of the sl{2) algebra is recovered by computing the corresponding Poisson 
bracket in (R 2 ,dgi A dpi). Moreover, its n-th coproduct yields a symplectic 
realization in R 2 ™ via 

^- = q 2 , Jo = P-q, J+ = p 2 + J2 b jiJ 2 > ( 7 ) 

3 

with bj 6 K. Again, the Lie-Poisson product is given by the standard Poisson 
bracket 

If \ = ^L^i_ _ df_dg_ 
<9q dp dp <9q 

and the Casimir ([5]) is simply 

6,q* 



where L 2 denotes the angular momentum. From this expression it is apparent 
that C is a homogenous function of degree 0. 
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This fit(2) ® ■ • • (gist(2) symmetry endows any ?i-dimensional Hamiltonian in 
the enveloping algebra with 2n — 3 integrals other than the Hamiltonian given 
by the left and right partial Casimirs. More precisely we have the following 
result, which we quote from Ref. [3]. 

Theorem 1. Any Hamiltonian Ho(p, q) = h{J + , J$,JJ) possesses 2n — 3 first 
integrals given by the left and right partial Casimirs 



E 

1<2<j <m 



(QiPj ~ QjPi) H — 



i=l 



C (m) - E 

n— m<i<j<n 



(qiPj - qjPi 



21 



E 



i— n— m+1 

where 1 < m < n. Here C^ n ' — Ci n ) = C is the Casimir ([6]) and the functions 

{H , C {l \C\ m) :l <l <n,l <m<n} 

are functionally independent. Moreover, the subsets {Ho, C^" 1 ^ : 1 < m < n} 
and {Ho, Cr m -\ : 1 < m < n} are in involution. 

The previous result obviously applies to the Hamiltonian which corre- 
sponds to 

TT = J +- C 

2(k + J_) ' 

The main result of this paper is that the remaining first integral, which makes 
H MS, can indeed be found easily. 



Theorem 2. The remaining first integral can be chosen as 

Ii(p, q) = Pi - 2-H"(p, q) qf + hq- 2 , 
for any 1 < i < n. Moreover, the set {Jj : 1 < i < n} is also in involution. 
Proof. The equations of motion under H are given by 



(8) 



Pi 



Pi 



2H qi + btq- 



k -\- q 

Combining both equations one can write 



(9) 
(10) 



O + q 2 )* +2(q- q)?» 



2Hq, + b t q: 



= 0. 



Multiplying this equation by 2(k + q 2 )<7i one immediately finds 
d 



dt 



[(K + q 2 ) 2 q 2 -2Hql + b iq r 2 ] =0, 



4 



which yields the desired result. The functional independence of Ii is easily 
established through a tedious but straightforward computation. The fact that 
{Ii} are in involution follows from 

{Ii,Ij} = -2q 2 ({H,p 2 } - 2H{H,q?} + b 3 {H, q r 2 }) 

+ 2q 2 ({H,p 2 } - 2H {H,q 2 } + h{H,qr 2 }) 

= -2q 2 {H,I ] } + 2q 2 {H 1 I l } = 0. 

□ 

Quadratic integrability is linked to the separability of the HJ equation by a 
theorem due to Kalnins and Miller [27]. In fact, it is not difficult to show that 
this equation is actually superseparable. If we write S(t, q) = W(q) — \Et the 
HJ equation for H reads 

3 

Thus one essentially recovers the HJ equation for the SW system with a different 
set of constants, namely 

#sw(p, q) - P 2 - Eq 2 + bjqj 2 ■ 

3 

Therefore, 

Proposition 3. H separates in the same coordinate systems as the SW system. 

In three dimensions, e.g., it separates in 8 out of the 11 possible orthogonal 
coordinate systems [16,23]. As a matter of fact, it has long been conjectured 
that any MS system should separate in multiple coordinate systems, but to our 
best knowledge this general claim has not been proved or disproved. 

3 Integration of the equations of motion 

In this section we shall compute the trajectories of the Hamiltonian system @ 
in closed form. It is convenient to start by exploiting the sl(2) symmetry to 
compute the evolution of the radial variable 

x = k + J_ = K + q 2 . 

We use the notation 

E:=2H= J -±^- 

x 

and assume c ^ 0. 

The evolution of the sl(2) generators (O is given by the following set of 
equations: 

x = {J-,H} 
Jo = {Jo, H} 
j+ = {J+,H} 



2J Q 



J+ + E(x - k) 



2EJ 



(11a) 
(lib) 
(11c) 
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Using these equations one can write the Casimir © as 

C = -\x 2 x 2 + Ex 2 + (c - Ek)x - ck > . 
For the sake of concreteness, let us assume that E > and set 

^■■=a( k+ e) + e- < 13 > 

In this case 

x 2 x 2 = AE [{x - a) 2 - 7 2 ] , (14) 
and one can easily integrate this equation as 

± 2VE(t -t) = ^{x - a) 2 - 7 2 + a cosh" 1 (^^) • (15) 

Here r is an arbitrary constant and 7 is the positive square root of 7 2 . Note 
that the equation (fT4|) for the variable x (the squared radius) coincides with the 
radial equation in the Kepler problem. 

The above expression yields t as a monotonic function of x, so this relation 
is globally invertible (in each half-orbit). However, it is not possible to obtain 
the inverse function x(t) in closed form. Hence we prefer to parametrize the 
trajectory by the radial variable x. In terms of this variable, the first integral ([5]) 
reads 

h = AE[(x a) 2 - 7 2 ] (g) 2 - Eq 2 + b iQ f . (16) 

Let us set 

Qi = qf, OLi = j 2 = a 2 + E~ l b h (17) 

so that 

Eqf + hq 2 -b i = E[(Qi - a,) 2 - jf] . 
Then one can perform the integration of Eq. (fTr?)) to obtain 

dQi f dx 

<pi ' 



, 2 



i.e., 

-1 f Qi - a 



cosh 



where y>j are constants and 7^ > 0. The full solution of this equation is therefore 
given by 



= on + 7i cosh(X + <pi) 

on + 7 _1 7i cosh(/3j (x - a) + 7* sinh^j 







1 - 


( £ T £ ) i 



(18) 

with X := cosh _1 (7~ 1 (a; — a)). Eqs. (fT5"]) and (jTS]) are the main results of this 
section. 
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All the orbits with positive energy are recovered through an appropriate 
choice of the In parameters 

{r, Oj, tpj : 1 < i < n - 1, 1 < j < n} , (19) 

where r was defined in Eq. (|15p . In what follows we shall express the remaining 
quantities in terms of this fundamental set. 

If we sum over i in Eq. (|18[) and use that x — n + ^\ Qi, we immediately 
find 

x — k = cti + "f^ 1 (x — a) ji cosht^i + 1— ( ^ 7f sinh ifj . 

i i i 

Therefore we reach the following compatibility conditions: 

on + k = a , (20a) 

i 

7i cosh ipi = j, (20b) 
y 7, sinh yjj = . (20c) 

Eq. (|20ap gives the value of a in terms of {ctj} and Eq. (|20b[) consistently 
provides the value of the Casimir by means Eq. ([T5)) . Eq. (|20c|) . which imposes 
the constraint 

^sinh^ (a? + S- 1 ^,;) 172 = 

on the the parameters {c^, by virtue of Eq. (|17|) . is responsible for the fact 
that only 2n — 1 among these parameters have been included in the set (|19p . 
Note that Eqs. (ITS]) and (|20a[) can be combined to express E as a function of 
{aj as 

Finally, this permits to express 7, in terms of {<x, } by means of Eq. (fTTj) . We 
omit the discussion of the case E < 0, which goes along the same lines. 



4 The geometric content of 7i 

It is apparent that TC is a (nonconstant) multiple of the SW system. What we 
want to stress in this section is that their connection does not end here: H. is in 
fact the intrinsic SW Hamiltonian in the manifold A4 n . And what makes this 
space remarkable is that its intrinsic SW model is also MS. 

It is apparent from Eq. that the Riemannian manifold A4 n is spherically 
symmetric, this SO(n) symmetry being a consequence of the sl(2) coalgebra 
structure outlined in Sec. [2] Its curvature is certainly nonconstant; for the sake 
of completeness we note that its scalar curvature is negative and given by 

1 , 3(n-2)q 2 + 2^n 

R = —{n— 1) -. ^5 . 

(K + q 2 ) 3 
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One can define intrinsic versions of the Kepler and harmonic oscillator po- 
tentials on A4 n . To this end, denote by A M 3 the Laplace-Beltrami operator 
on M 3 and let V3 € C°°(A1 3 \{0}) be its (minimal symmetric) Green function, 
i.e., a function V^q) = v(\q\) satisfying 



-& M3 v 3 



Here v £ C°°(]R + \{0}) and Sj^s stands for the delta distribution in Ai 3 sup- 
ported at the origin. It can be proved [13] that such a function exists and is 
unique. 

Definition 4. The Kepler and harmonic oscillator potentials in M. n are 

^Kepler(q) := A"«(|q|) , V Harm (q) := K v(\q\)- 2 , 
where K € R is an arbitrary constant. 

Remark 5. The above definition, which is based on the potential theory of these 
spaces, can be easily extended to any spherically symmetric Riemannian space 
with the appropriate behavior at infinity [30] . Such a prescription reproduces the 
intrinsic Kepler and harmonic potentials on constant curvature spaces (see [12] 
and references therein). The case of U(2)-symmetric Kahler 4-manifolds, which 
had already appeared in the literature [33] , is based on the same ideas but does 
not fit into this framework. Certainly the above prescription does not generally 
lead to MS potentials. 

Let q € M 3 and define r — |q|. It can be readily verified that the action of 
the Laplace-Beltrami operator in M 3 on some function f{r) is 

df 



M ^ r 2 (n + r 2 )dr\ dr 

By setting the above expression equal to zero for r > it is straightforward to 
find that the function v defined above is 



\Jk + r 2 
v[r) = 



up to a multiplicative constant. Hence the Kepler and harmonic oscillator po- 
tentials in M n are respectively given by 



XT , X „ V K + q T 7 , S Kq 2 

VKepler(qJ = K V , Vkarm(q) = ; o • 

|q| K + q 2 

This shows that the MS Hamiltonian ([T|) is in fact the intrinsic SW systems, 
i.e., 

n(p, q) = ||p||^„ + Vkarm(q) + (« + q 2 )" 1 h ilf > 

3 

with || • \\m™ representing the metric on the cotangent bundle of M n . 



8 



Acknowledgements 



This work was partially supported by the Spanish MEC and the Junta de 
Castilla y Leon under grants no. FIS2004-07913 and VA013C05 (A.B. and 
F.J.H.), by the Spanish DGI under grant no. FIS2005-00752 (A.E.) and by 
the INFN-CICyT (O.R.). Furthermore, A.E. acknowledges the financial sup- 
port of the Spanish MEC through an FPU scholarship, as well as the hospitality 
and the partial support of the Physics Department of Roma Tre University. 

References 

[1] M. A. Agrotis, P. A. Damianou, C. Sophocleous, The Toda lattice is super- 
intcgrablc, Physica A 365 (2006) 235-243. 

[2] V. I. Arnold, V. V. Kozlov, A. I. Ncishtadt, Mathematical aspects of clas- 
sical and celestial mechanics, Springer- Verlag, Berlin, 1997. 

[3] A. Ballesteros, F. J. Herranz, Universal integrals for supcrintcgrable sys- 
tems on A-dimensional spaces of constant curvature, J. Phys. A: Math. 
Theor. 40 (2007) F51-F59. 

[4] A. Ballesteros, F. J. Herranz, M. Santander, T. Sanz-Gil, Maximal super- 
intcgrability on TV-dimensional curved spaces, J. Phys. A: Math. Gen. 36 
(2003) L93-L99. 

[5] A. Ballesteros, F. Musso, O. Ragnisco, Classical and quantum integrable 
systems: the coalgebra approach, Rcgul. Chaotic Dyn. 4 (2002) 393-398. 

[6] A. Ballesteros, O. Ragnisco, A systematic construction of completely in- 
tegrable Hamiltonians from coalgebras, J. Phys. A: Math. Gen. 31 (1998) 
3791-3813. 

[7] J. Bertrand, Mecanique analytique, C. R. Acad. Sci. 77 (1873) 849-853. 

[8] A. T. Bruce, R. G. McLenaghan, R. G. Smirnov, A geometrical approach 
to the problem of integrability of Hamiltonian systems by separation of 
variables, J. Geom. Phys. 39 (2001) 301-322. 

[9] F. Calogcro, Solution of the one-dimensional A-body problems with 
quadratic and/or inversely quadratic pair potentials, J. Math. Phys. 12 
(1971) 419-436. 

[10] F. Calogero, Two new classes of isochronous Hamiltonian systems, J. Non- 
lin. Math. Phys. 11 (2004) 208-222. 

[11] J. F. Carinena, M. F. Ranada, M. Santander, A nonlinear deformation of 
the isotonic oscillator and the Smorodinski-Winternitz system: integrabil- 
ity and superintegrability, Regul. Chaotic Dyn. 10 (2005) 423-436. 

[12] J. F. Carinena, M. F. Ranada, M. Santander, Central potentials on spaces 
of constant curvature: the Kepler problem on the two-dimensional sphere 
S 2 and the hyperbolic plane H 2 , J. Math. Phys. 46 (2005) 52702(26). 



9 



[13] A. Enciso, D. Peralta-Salas, Geometrical and topological aspects of elec- 
trostatics on Riemannian manifolds, submitted for publication. 

[14] N. W. Evans, Superintegrability in classical mechanics, Phys. Rev. A 41 
(1990) 5666-5576. 

[15] N. W. Evans, Superintegrability of the winternitz system, Phys. Lett. A 
147 (1990) 483-486. 

[16] N. W. Evans, Group theory of the Smorodinsky-Winternitz system, J. 
Math. Phys. 32 (1991) 3369-3375. 

[17] I. Fris, V. Mandrosov, J. Smorodinsky, M. Uhlir, P. Winternitz, On higher- 
order symmetries in quantum mechanics, Phys. Lett. 16 (1965) 354-356. 

[18] C. Gonera, Isochronic potentials and new family of superintegrable systems, 
J. Phys. A: Math. Gen. 37 (2004) 4085-4095. ' 

[19] Y. I. Granovski, I. M. Lutzenko, A. S. Zhedanov, Mutual integrability, 
quadratic algebras, and dynamical symmetry, Ann. Phys. 217 (1992) 1-20. 

[20] S. Gravel, P. Winternitz, Superintegrability with third-order integrals in 
quantum and classical mechanics, J. Math. Phys. 43 (2002) 5902-5912. 

[21] F. J. Herranz, A. Ballesteros, Superintegrability on three-dimensional Rie- 
mannian and relativistic spaces of constant curvature, SIGMA 2 (2006) 
010(22). 

[22] A. Huaux, Sur la separation des variables dans l'equation aux derivees 
partielles de Hamilton-Jacobi, Ann. Mat. Pura Appl. 108 (1976) 251-282. 

[23] E. G. Kalnins, Separation of variables for Riemannian spaces of constant 
curvature, Wiley, New York, 1986. 

[24] E. G. Kalnins, J. M. Kress, W. Miller, Second-order superintegrable sys- 
tems in conformally flat spaces. IV. The classical 3D Stackel transform and 
3D classification theory, J. Math. Phys. 44 (2006) 43514(40). 

[25] E. G. Kalnins, J. M. Kress, W. Miller, P. Winternitz, Superintegrability 
in a two-dimensional space of nonconstant curvature, J. Math. Phys. 43 
(2002) 970-983. 

[26] E. G. Kalnins, J. M. Kress, W. Miller, P. Winternitz, Superintegrable sys- 
tems in Darboux spaces, J. Math. Phys. 44 (2003) 5811-5848. 

[27] E. G. Kalnins, W. Miller, The theory of orthogonal i?-separation for 
Helmholtz equations, Adv. Math. 51 (1984) 91-106. 

[28] E. G. Kalnins, G. C. Williams, W. Miller, G. S. Pogosyan, On superinte- 
grable symmetry-breaking potentials in iV-dimensional Euclidean space, J. 
Phys. A: Math. Gen. 35 (2002) 4755-4773. 

[29] G. Kocnigs, Sur les geodesiques a integrates quadratiques, in: G. Darboux 
(Ed.), Lecons sur la theorie generale des surfaces, Vol. 4, Chelsea, New 
York, 1972, pp. 368-404. 



10 



[30] P. Li, L. F. Tarn, Symmetric Green's functions on complete manifolds, 
Amer. J. Math. 109 (1987) 1129-1154. 

[31] H. Liebmann, Nichteuklidische Geometrie, Walter de Gruyter, Berlin, 1923. 

[32] J. Moser, Three integrablc Hamiltonian systems connected with isospectral 
deformations, Adv. Math. 16 (1975) 197-220. 

[33] A. Nersessian, A. Yeranyan, Three-dimensional oscillator and Coulomb sys- 
tems reduced from Kahler spaces, J. Phys. A: Math. Gen. 37 (2004) 2791- 
2801. 

[34] M. F. Ranada, Superintegrable n = 2 systems, quadratic constants of mo- 
tion and the potentials of Drach, J. Math. Phys. 38 (1997) 4165-4178. 

[35] M. B. Sheftel, P. Tcmpcsta, P. Winternitz, Superintegrable systems in 
quantum mechanics and classical Lie theory, J. Math. Phys. 42 (2001) 659- 
673. 

[36] B. Sutherland, Exact results for a quantum many-body problem in one 
dimension, Phys. Rev. A 4 (1971) 2019-2021. 

[37] M. Toda, One-dimensional dual transformation, J. Phys. Soc. Japan 22 
(1967) 431-436. 

[38] S. Wojciechowski, Superintegrability of the Calogero-Moser system, Phys. 
Lett. A 95 (1983) 279-281. 



11 



